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ABSTRACT 

An abundance analysis is presented for IRAS 18095+2704 (V887 Her), a post-AGB star 
and proto-planetary nebula. The analysis is based on high-resolution optical spectra from the 
McDonald Observatory and the Special Astrophysical Observatory. Standard analysis using a 
classical Kurucz model atmosphere and the line analysis program MOOG provides the atmo- 
spheric parameters: T e g = 6500 K, log g = +0.5, and a microturbulent velocity £ = 4.7 km 
s _1 and [Fe/H] = —0.9. Extraction of these parameters is based on excitation of Fel lines, 
ionization equilibrium between neutral and ions of Mg, Ca, Ti, Cr, and Fe, and the wings 
of hydrogen Paschen lines. Elemental abundances are obtained for 22 elements and upper 
limits for an additional four elements. These results show that the star's atmosphere has not 
experienced a significant number of C- and s-process enriching thermal pulses. Abundance 
anomalies as judged relative to the compositions of unevolved and less-evolved normal stars 
of a similar metallicity include Al, Y, and Zr deficiencies with respect to Fe of about 0.5 dex. 
Judged by composition, the star resembles a RV Tauri variable that has been mildly affected 
by dust-gas separation reducing the abundances of the elements of highest condensation tem- 
perature. This separation may occur in the stellar wind. There are indications that the standard 
ID LTE analysis is not entirely appropriate for IRAS 18095+2704. These include a supersonic 
macroturbulent velocity of 23 km s _1 , emission in Ha and the failure of predicted profiles to 
fit observed profiles of H/3 and H7. 
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1 INTRODUCTION 

Stars of low and intermediate mass with initial masses between 
O.&Mq - 8M0 evolve to the asymptotic giant branch (AGB). Then, 
thanks to severe mass loss, the AGB star evolves rapidly at nearly 
constant luminosity to higher effective temperatures to the white 
dwarf cooling track. Typical stellar lifetimes of post-AGB stars are 
expected to be of the order of 10 4 years (Schonberner 1983). The 
gas lost by the AGB star forms a circumstellar shell. When the post- 
AGB star is cool, the dust in the shell heated by stellar radiation 
provides an infrared excess. When the star has traversed the top of 
the H-R diagram to higher effective temperatures, the circumstellar 
gas is ionized. Then, the star is said to have evolved to reached the 
proto-planetary nebula stage. Shortly after this, the post-AGB star 
has evolved to become a planetary nebula with a hot white dwarf 
as the central star. Determinations of the chemical composition for 
post-AGB star hold the potential of yielding insights into the chem- 
ical history of the AGB star and its conversion by mass loss to its 
slimmer post-AGB form. 
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In this paper, we present a determination of the chemical compo- 
sition of IRAS 18095+2704, a post-AGB star with a substantial 
dusty circumstellar shell. The discovery of the optical counterpart 
IRAS 18095+2704 was made by Hrivnak, Kwok, & Volk (1987, 
1988). This V = 10.4 mag star is a high-latitude F supergiant with 
a large far-IR excess. In the Catalog of Low Resolution IRAS Spec- 
tra assembled by Hrivnak, Kwok, & Volk (1988), the star has a 
peculiar IR continuum slope at wavelengths shortward of the 10 
(jrn silicate emission feature. According to Volk & Kwok (1987), 
this peculiar continuum shape is a result of a detached dust shell. 
Observational evidence for an expanding shell came from Lewis, 
Eder, & Terzian (1985) and Eder, Lewis, & Terzian (1988) via de- 
tection of OH maser emission at 1612 and 1665/67 MHz from the 
Arecibo telescope. Gledhill et al. (2001) from imaging polarimetry 
report an extended envelope or a reflection nebula around the star. 



The pioneering study of IRAS 18095+2704's composition was re- 
ported by Klochkova (1995) from echelle spectra (R = 24 000) 
covering the wavelength ranges 5050 A to 7200 A and 5550 A to 
8700 A. Abundances of 26 elements were obtained. The star was 
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Figure 1. The spectrum for IRAS 18095+2704 over the wavelength region 4513-4592 A . The McDonald spectrum is plotted in black and the SAO spectrum 
in green. Selected lines are identified. 



found to be moderately metal poor, [Fe/H]= —0.780 with a relative 
enrichment of C and N, i.e., [C/Fe]=0.5 and [N/Fe]=0.5, as might 
be expected of a post-AGB star evolved from a C-rich AGB star. 
Although other elements up through the Fe-peak had roughly their 
anticipated abundances, two results drew comment. First, there was 
a difference in abundances [X/H] derived from neutral and first- 
ionized lines of several elements, i.e., differences of 1.6 (Ti), 1.4 
(V), 1.0 (Cr), and 2.2 dex (Y). For Fe, this difference was zero 
because it was the condition enforced in determining the surface 
gravity. Second, the accessible lanthanides (La, Pr, Nd and Eu) 
represented by ionized lines were overabundant by about [X/Fe] 
~ +0.7 relative to what is expected for an unevolved metal-poor 
star. Although one might attribute this overabundance to s-process 
enrichment expected of a post-AGB star, one notes that Eu, pre- 
dominantly an r-process element, had the highest overabundance 
([Eu/Fe] = +1.2), Ba was not overabundant ([Ba/Fe] = —0.2), and 
yttrium, the sole representative of lighter s-process species, as ana- 
lyzed from Y II lines, was also not overabundant ([Y/Fe] ~ 0). 

In this paper, we determine afresh the composition of 
IRAS 1 8095+2704 from echelle optical spectra: one spec- 
trum was obtained at the W.I. McDonald Observatory and another 
at the Special Astrophysical Observatory (SAO). In addition to ex- 
amining the unusual results noted above, we seek an interpretation 
of the star's the composition in light of its proposed status as a 
slightly metal-poor post-AGB star. 



250 



§ 200 
o 

n 



< 150 

EC 




50 100 150 200 250 

EW(mA)(McD) 



Figure 2. Equivalent widths (EW) for 69 lines in the McD spectrum versus 
the EWs of the same lines from SAO (2009) spectrum. The solid line is the 
least-squares fit. The dash-dot line indicates EW quality. 



1 Standard notation is used for quantities [X] where [X]=log(X) s t a 
log(X) . 



2 OBSERVATIONS AND DATA REDUCTION 

The McDonald spectrum for the abundance analysis was obtained 
on the night of 2008 August 10 (JD 2454688.7) at the 2.7 me- 
ter Harlan J. Smith reflector with the Tull cross-dispersed echelle 
spectrograph (Tull et al. 1995) at a spectral resolution of X/dX ~ 
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Table 1. Absorption components of Nal Di (5895.923 A ) and Nal D2 
(5889.953 A ) lines in the McD spectrum of IRAS 18095+2704 (V887 Her). 
Vq are the heliocentric radial velocities of the components. 





NaDl 


NaD2 


Component 


Vq 


Vq 




(kms 1 ) 


(kms 1 ) 


1 


-43.6 


-41.7 


2 


-27.1 


-26.1 


3 


-18.1 


-17.0 


4 


-6.8 


-4.5 



60, 000. The spectrum covers the wavelength ranges 3800 A to 
10500 A with no gaps in the wavelength ranges 3800 A to 4885 
A and 5020 A to 5685 A but coverage is incomplete but substantial 
beyond 5700 A; the effective short and long wavelength limits are 
set by the useful S/N ratio. A ThAr hollow cathode lamp provided 
the wavelength calibration. Flat-field and bias exposures completed 
the calibration files. The signal-to-noise ratio ranges between 40 
and 90 per pixel, not only changing with the blaze function within 
echelle orders but also star brightness between echelle orders. 
The McD observations were reduced using the STARLINK 
echelle reduction package ECHOMOP (Mills & Webb 1994). 
The spectra were extracted using ECHOMOP' S implementation 
of the optimal extraction algorithm developed by Home (1986). 
ECHOMOP propagates error information based on photon statis- 
tics and readout noise throughout the extraction process. The bias 
level in the overscan area was modeled with a polynomial and sub- 
tracted. The scattered light was modeled and removed from the 
spectrum. In order to correct for pixel-to-pixel sensitivity varia- 
tions, 'flatfield' exposures from a halogen lamp were used. Indi- 
vidual orders were cosmic-ray cleaned, and continuum normalized 
with bespoke echelle reduction software in IDL (§ahin 2008). Re- 
duced spectra were transferred to the STARLINK spectrum analy- 
sis program DIPSO (Howarth et al. 1998) for further analysis (e.g. 
for equivalent width measurement). In equivalent width measure- 
ments, local continua on both side of the lines were fitted with 
a first-degree polynomial then equivalent widths were measured 
with respect to these local continua using a fitted Gaussian profile. 
For strong lines, a direct integration was preferred to the Gaussian 
approximation. The errors for each equivalent width measurement 
were determined on the basis of scatter of linear continuum fit and 
signal-to-noise ratio of each measured line in the spectra. Errors on 
the measured equivalent widths are calculated using the prescrip- 
tions given by Howarth & Phillips (1986). 

The SAO spectrum was obtained on the night of 2009 June 10 (JD 
2452993.4) by VK and NST with the NES echelle spectrograph 
mounted at the Nasmyth focus of the 6-m telescope of the Special 
Astrophysical Observatory (Panchuk et al. 2007) with a 2048 X 
2048 CCD with an image sheer (Panchuk et al. 2007) and a spec- 
tral resolution of X/d\ ^ 60,000. A modified ECHELLE con- 
text (Yushkin & Klochkova 2005) of MIDAS package was used to 
extract one-dimensional vectors from the two-dimensional echelle 
spectra. Wavelength calibration was performed using a hollow- 
cathode Th-Ar lamp. The wavelength coverage for the SAO spec- 
trum was 4460 - 5920 A . Measurement of equivalent widths was 
carried out as for the McDonald spectrum. 

The agreement between the McDonald and SAO spectra is sat- 
isfactory, i.e., the two spectra are similar as to line width, depth and 
equivalent width for weak to strong lines. This is shown by the 



section of the reduced spectra illustrated in Figure [T] A sample of 
apparently unblended lines was selected from across the common 
wavelength interval and their equivalent widths (EWs) measured in 
both the McDonald and the SAO spectra. The comparison of EWs 
shown in Figure [2] shows good agreement between the two sets of 
measurements. Across the common wavelength interval, we com- 
pare EWs from McDonald and SAO spectra, especially for lines at 
the limit of detection and for elements represented by just one or 
two lines. The SAO spectrum was used to provide lines that fell in 
the inter-order gaps of the McDonald spectrum. 



3 GENERAL FEATURES OF THE SPECTRA 

Although the data are sparse, the star is probably not a large am- 
plitude velocity variable. Klochkova (1995) reported a heliocen- 
tric radial velocity of — 32.5±0.4 km s - . Hrivnak, Kwok, & Volk 
(1988) report — 30±2 km s~ from an unspecified number of mea- 
surements but add that there is 'an indication of variability'. The 
McD spectrum gives — 30±1 km s _1 from the metal lines. The 
2009 SAO spectrum gives -31.8±1.7 km s _1 . Although Lewis, 
Eder, & Terzian (1985) and Eder, Lewis, & Terzian (1988) cite the 
heliocentric velocity as —17.4 km s _1 from their observed OH ra- 
dio lines, Hrivnak (2009, private communication) indicates that this 
value resulted from an incorrect conversion of LSR to heliocentric 
velocity and a velocity of about —30 km s _1 is obtained from the 
OH velocities. 

The Na D lines show four components. Figure [3] shows the Na D 
from the McDonald and SAO spectra. Heliocentric velocities of the 
four principal components in the McDonald spectrum are listed in 
Table 1. The SAO spectrum gives similar velocities. Stellar photo- 
spheric Na lines at about —30 km s~ x must be largely masked by 
these multiple circumstellar components. Component 1, if not an 
interstellar component, represents an outflow at a velocity of about 
12 km s . Components 3 and 4 are falling toward the star at ve- 
locities of about 12 and 25 km s _1 , respectively. 

Weak emission in the blue and red wings of the Ha profile flanking 
a deep narrow absorption core was reported by Klochkova (1995 
- see also Tamura, Takeuti & Zalewski 1993). On the McDonald 
spectrum, Ha occurs at the very edge of an order but a very shallow 
a deep core flanked by red emission is seen. H/3 and higher lines in 
the Balmer series and Paschen lines are purely in absorption. The 
2009 SAO spectrum did not include Ha. 

The stellar absorption lines are broad. If one accepts classical no- 
tions of microturbulence and macroturbulence, this width suggests 
substantial macroturbulence in the atmosphere. The microturbu- 
lence is about 5 km s _1 (see below). The instrumental width is 
about 5 km s~\ Correcting for the microturbulence and the instru- 
mental width, the macroturbulence is estimated to be about 23 km 
s _1 . This represents a highly supersonic velocity. 

Strong lines clearly show asymmetric profiles with an extended 
blue wing. This is illustrated in Figure [4] where the O I triplet lines 
at 777 1-7775 A and the Sill 6347A line are shown. The blue wing 
extension extends from about —25 km s _1 to —50 km s _1 with re- 
spect to the photospheric velocity. Inspection of strong unblended 
lines shows that the blue asymmetry is present also for low exci- 
tation strong lines (e.g., Mglb 5167 A , 5172 A , and SrII 4215 
A ). This outward motion may represent the early stages of a stel- 
lar wind or inhomogeneities (i.e., stellar super-granulation) in the 
photosphere. 
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Figure 3. The profile of the D-lines of Nal a) in the McDonald 2008 spectrum and b) in the SAO 2009 spectrum of IRAS 18095+2704. Telluric H 2 lines 
have not been removed. 
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Figure 4. Profiles of Si II 6347 A line and O I 7772 A triplet in the McDon- 
ald spectrum. The best-fit synthetic spectra (in green) fail to account for the 
extended blue wings. The bottom of each panel shows the residuals shifted 
by 0.3 relative flux for clarity. 
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Figure 5. Determination of atmospheric parameters T e g and £ using abun- 
dance (log e(Fe)) as a function of both lower level excitation potential 
(LEP) and reduced equivalent width (log(W/ A)). In all the panels, the solid 
black line is the least-square fit to the data. The zero slope of this line is for 
T c ff =6500 K. The ionized (in red) iron lines were over plotted. 



4 ABUNDANCE ANALYSIS - THE MODEL 

ATMOSPHERES AND STELLAR PARAMETERS 

The abundance analysis was performed using the local thermo- 
dynamic equilibrium (LTE) stellar line analysis program MOOG 
(Sneden 2002). Model atmospheres were obtained by interpolating 
in the ATLAS9 model atmosphere grid (Kurucz 1993). The models 
are line-blanketed plane-parallel uniform atmospheres in LTE and 



hydrostatic equilibrium with flux (radiative plus convective) con- 
servation. A model is defined by an effective temperature T c ff, sur- 
face gravity g, chemical composition as represented by metallicity 
[Fe/H] and a microturbulence velocity £ of 2 km s _1 . 

Several of the assumptions adopted in the construction and appli- 
cation of the model atmospheres are of uncertain validity when 
considering post-AGB stars and supergiants in general. The pres- 
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Figure 6. The standard deviation of the Fe and Ti abundances from the suite 
of Fe I and Ti II lines as a function of the microturbulence £. 
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Figure 7. The T c ff versus log g plane showing the various loci discussed in 
the text. Loci include those derived from the Paschen line profiles (Pn,i4, 
dashed line), the excitation of the Fe I lines (Fe I(ext), dashed line), and the 
imposition of ionization equilibrium for Mg, Si, Ca, Cr, and Fe. 



ence of supersonic levels of macroturbulence and hints of a stellar 
wind seem incompatible with the assumption of hydrostatic equi- 
librium. Departures from LTE are probable for these low density 
atmospheres. The real atmosphere is likely to depart from the uni- 
form plane-parallel layered theoretical construction. These qualify- 
ing remarks should be borne in mind when interpreting the derived 
abundances. 



4.1 Spectroscopy - Fe I and Fe II lines 

A standard spectroscopic method of determining the effective tem- 
perature, surface gravity, and the microturbulence using Fel and 
Fe II was applied. Application of strict criteri43 in the selection of 
suitable lines provided a list of 31 Fe I lines with lower excitation 
potentials (LEP) ranging from 0.9 to 4.5 eV and EWs of up to 224 
mA and 1 1 Fe II lines with excitation potentials of 2.8 to 3.9 eV and 
EWs of up to 154 mA. Only six Fe I lines have EW greater than 150 
mA. These lines are listed in Table 2. The ^/-values are taken from 
Fiihr & Wiese (2006). 

One estimate of the temperature is found from Fe I lines by exci- 
tation balance. The value for T e g is chosen so that the abundance 
is independent of a line's lower level excitation potentials (LEP). 
The microturbulence is determined by adjusting the £ so that the 
abundances were independent of reduced equivalent width (W/X). 
For our sample of Fe I lines, these two conditions are imposed si- 
multaneously (see Figure 5). The microturbulence may also be de- 
termined from the Ti II lines as they are numerous and of suitable 
strength. For a given model, we compute the dispersion in the Fe 
(or Ti) abundances over a range in the £ from 2 to 10 km s _1 . In 
Figure|6] the dispersion a for Fe and Ti lines is displayed. The dis- 
persion a values for the Fe I and Ti II lines are computed for two 
different effective temperature values: T e « = 5750 K (e.g. green 
curve for Fe) and 6500 K (e.g. blue curve for Ti). From the dis- 
persion g vs microturbulence velocity £ plots for the Fe lines, the 
microturbulence velocity is found to be in the range 4.7^^6.0 km 
s _1 . A minimum value of a for the Ti lines is reached at £ ~ 5.0 
km s _1 . We adopt 4.7±0.5 km s _1 . The solutions for £ are not 
particularly dependent on the chosen T e g of the model. 

Imposition of ionization equilibrium for an element represented by 
lines from neutral atoms and singley-charged ions provides a locus 
in the temperature-gravity plane running from low T e g and low g to 
high T ff and high g. Such loci are shown in Figure|7]for Mg, Si, Ca, 
Cr, and Fe; silicon presents a problem (see below). These loci with 
the (vertical) loci provided by the T e g from excitation of Fe I lines 
and the Paschen line profiles (see below) give T e g =6500 K and 
log g =+0.5 cgs when greatest weight is given to the Fe ionization 
equilibrium on account of the greater number of Fe lines. 



4.2 The Balmer and Paschen lines 

In principle, Balmer line profiles offer an alternative method of es- 
timating atmospheric parameters. For warm supergiants like IRAS 
18095+2704, the Balmer lines are sensitive to T e g and insensitive 
to log p so that they provide an isothermal in Figure [7] Ha is un- 
suited to this purpose because it shows strong emission distorting 
the expected deep and broad photospheric absorption profile. Both 
H/3 and H7 appear to have the anticipated profiles. Predicted pro- 
files are computed with SYNTHE as MOOG does not compute 
synthetic profiles for the hydrogen lines. Figure [8] shows observed 
and predicted profiles from synthetic spectrum calculations for T c s 
from 5500 K to 6500 K. It is apparent that for both lines the best- 
fitting predicted profile corresponds to an effective temperature of 
around 5750 K; this estimate is almost independent of the adopted 
surface gravity. This temperature is significantly cooler than that 
suggested by the intersection of the Fe excitation temperature and 
the various ionization equilibria (Figure|7]). 



Lines that are well defined and not blended. 
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Figure 8. Observed and model line profiles for H/3 and H7. The green, red, blue, and black lines show the theoretical profiles for T c ff = 5500 K, 5750 K, 6000 
K, and 6500 K, respectively, all for a surface gravity log g = 0.5 dex. 
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Figure 9. Observed and model line profiles for hydrogen Paschen lines. The green and red lines show the theoretical profiles for T e ff = 5750 K, and 6500 K, 
respectively, both for a surface gravity log g = 0.5 dex. 
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Emission in the Balmer lines may be responsible in part for a 
systematically lower estimate of the effective temperature. The 
Ha profile is very strongly affected by emission at the time of 
our observation. In order to provide a more complete assessment 
of the emission, we use the Ha profile illustrated by Khochkova 
(1995). The observed profile is much shallower than a profile syn- 
thesized according to the adopted atmospheric parameters (i.e., 
left = 6500K). Obviously, emission is present across the entire 
profile and sharp emission just to the red of center hints at presence 
of a P Cygni profile. An emission profile expressed as the differ- 
ence between the adopted observed profile and that computed for 
6500 K has the shape necessary to correct the observed H/3 and H7 
profiles to the predicted profiles also for 6500 K. 

The outstanding issue is the question of the strength of the emis- 
sion at H/3 and H7. Standard Case A (where HI in the ionized gas 
is optically thin; i.e., of low density that Lyman photons can es- 
cape) and Case B (where HI in the ionized gas is optically thick to 
Lyman line photons, i.e., Lyman line photons in consequence get 
absorbed) decrements of the recombination theory for ionized gas 
suggest that flux ratios might be Ha/H/3 ~ 3 and H/3/Hy ~ 2.5 
(Osterbrock & Ferland 2006). Reduction of the Ha emission by 
these factors provides an emission component at H/3 and H7 that is 
slightly too weak to reconcile the observed (corrected for emission) 
with the profile predicted for 6500 K. Additionally, the fact that H7 
and H/3 require the same low effective temperature (T e fi — 5750 
K) suggests that superposition of emission on a photospheric profile 
may not be the explanation. Possibly, the Balmer line profiles are 
influenced by an upper photosphere (a highly-structured chromo- 
sphere with a developing wind?) that distorts the observed profiles 
in a way not modellable using Case A or Case B decrements. 

Paschen lines are expected to be less sensitive to structural adjust- 
ments of the upper photosphere and to emission. Figure [9] shows 
observed and synthesized profiles for a selection of Paschen lines 
recorded on the McDonald spectrum. Profiles for the model with 
T c ff = 6500 K, log p = +0.5 offer a fine fit to the observed pro- 
files. The latter are very poorly fit with the T e g — 5750K model 
suggested by the Balmer H/3 and H7 profiles. Thus, the Paschen 
profiles confirm the effective temperature from the excitation of 
the Fe I lines. The new isothermal provided by the hydrogen lines 
of the Paschen series at 8862 A (Pu) and 8598 A (P14) coincides 
with the locus from the excitation of Fe I lines in Figure UJ 
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Figure 10. The observed spectrum near 4774 A and synthetic spectra for 
the three carbon abundances shown on the figure. 
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5 ABUNDANCE ANALYSIS - ELEMENTS AND LINES 

For the lines of neutral and/or singly-ionized atoms, we conducted a 
systematic search by using lower excitation potential and gf -value 
as the guides. The REVISED MULTIPLET TABLE (RMT) (Moore 
1945) is used as an initial guide in this basic step. Lines chosen by 
this search are listed in Table 3. When a reference to solar abun- 
dances is necessary in order to convert our abundance of element X 
to either of the quantities [X/H] or [X/Fe], Asplund et al. (2009) is 
preferred. 

Comments on individual elements with notes on the adopted gf- 
values follow: 

C: The ^/-values are taken from Wiese, Ftihr & Deter (1996). Ten 
lines give a mean abundance log e(C)=7.92±0.17. Observed and 
synthetic spectra for three different carbon abundances are shown 
in Figure[l0]for a region providing three of the ten chosen lines. 



Figure 11. The observed spectrum near 8216A and synthetic spectra for the 
three nitrogen abundances shown on the figure. 

N: The (//-values are taken from Wiese, Fiihr & Deter (1996). 
Promising Nl lines are in the red outside the wavelength range 
covered by the SAO spectrum. The McDonald spectrum covers the 
region spanned by several multiplets. Several potential lines fall in 
inter-order gaps. Figure QT] shows a possible detection of one N I 
line. The N abundance is log e(N)^ 7.1 but this might properly be 
considered an upper limit. 

O: The (//-values for the permitted and forbidden Ol lines are 
taken from Wiese, Fiihr & Deter (1996). On the McDonald spec- 
trum, the forbidden oxygen lines at 5577 A, 6300 A and 6363 A 
are detected. Weak permitted lines of RMT 10 near 6156A are also 
analyzed. Figure [121 shows the best-fitting synthetic spectrum for 
the Ol 6156 A region. Forbidden and permitted lines give a very 
similar abundance. The strong O I triplet at 7774 A and the 8446 A 
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Table 2. Fe I and Fe II lines used in the analysis and corresponding abundances for a model of T c ff = 6500 K, log g = 0.5. RMT is the multiplet number 
given by the Revised Multiplet Table (Moore 1945). 



Species 


X 


EW loge(Fe) LEP 


log(9/) 


RMT 


Species 


A 


EW lo 


ge(Fe) LEP 




RMT 




(A) 






v ev ) 










(A) 


(mA) 


(dex) 


(eV) 






Fe I 


4447.717 


LOO 


0. / 4 


z.zz 


— 1 


Q/l O 

o*±A 


68 


Fel 


5429.697 


220 


P. BQ 
O.OO 





96 


— 1.879 


15 


Fel 


4459.118 


LOO 


P. TQ 
0. 1 O 


Z. lo 


1 

— 1 


07(1 

a i y 


68 


Fel 


5434.524 


111 


P A O 

0.4Z 


1 


01 


—2.122 


15 


Fel 


4476.019 




0.00 


Z.oD 


n 
— li 




350 


Fel 


5446.917 


206 


(\ 7/1 
D. / 4 





99 


— 1.914 


15 


Fel 


4531.148 


i on 
LAV 


O.o4 


i /in 
1.4y 


— A 


loo 


39 


Fel 


5569.618 


79 


P. A Q 


3 


12 


—0.486 


686 


Fel 


4871.318 


Lifo 


O.ol 


Z.ol 


n 


ODA 


318 


Fel 


5572.842 


113 


P. A 7 
0.4 I 


3 


10 


—0.275 


686 


Fel 


4872.138 


139 


6.54 


2.88 


-0 


567 


318 


Fel 


5615.644 


167 


6.54 


3 


33 


0.050 


686 


Fel 


4903.308 


85 


6.46 


2.88 


-0 


926 


318 


Fel 


5624.542 


61 


6.53 


3 


42 


-0.755 


686 


Fel 


4918.993 


168 


6.52 


2.87 


-0 


342 


318 


Fel 


6230.710 


68 


6.36 


2 


56 


-1.281 


207 


Fel 


5014.941 


60 


6.53 


3.94 


-0 


303 


965 


Fel 


6393.601 


56 


6.43 


2 


43 


-1.580 


168 


Fel 


5074.748 


75 


6.84 


4.22 


-0 


230 


1094 


Fel 


6411.649 


53 


6.61 


3 


65 


-0.718 


816 


Fel 


5090.774 


38 


6.69 


4.26 


-0 


440 


1090 


Fell 


4472.929 


154 


6.74 


2 


84 


-3.530 


37 


Fel 


5216.274 


98 


6.66 


1.61 


-2 


150 


36 


Fell 


4893.820 


58 


6.64 


2 


83 


-4.270 


36 


Fel 


5232.940 


224 


6.86 


2.94 


-0 


057 


383 


Fe II 


5120.352 


63 


6.63 


2 


83 


-4.214 


35 


Fel 


5302.302 


72 


6.49 


3.28 


-0 


720 


553 


Fell 


5132.669 


70 


6.55 


2 


81 


-4.090 


35 


Fel 


5324.179 


148 


6.42 


3.21 


-0 


103 


553 


Fell 


5254.929 


143 


6.64 


3 


23 


-3.227 


49 


Fel 


5364.871 


77 


6.59 


4.45 





228 


1146 


Fell 


5256.938 


71 


6.79 


2 


89 


-4.250 


41 


Fel 


5367.467 


93 


6.49 


4.42 





443 


1146 


Fe II 


5325.553 


128 


6.50 


3 


22 


-3.220 


49 


Fel 


5369.962 


105 


6.45 


4.37 





536 


1146 


Fe II 


5432.967 


93 


6.68 


3 


27 


-3.629 


55 


Fel 


5383.369 


113 


6.35 


4.31 





645 


1146 


Fell 


6147.741 


104 


6.40 


3 


89 


-2.721 


74 


Fel 


5393.168 


67 


6.40 


3.24 


-0 


715 


553 


Fell 


6149.258 


106 


6.54 


3 


89 


-2.840 


74 


Fel 


5397.128 


192 


6.64 


0.92 


-1 


993 


15 


Fe II 


6416.920 


93 


6.47 


3 


89 


-2.880 


74 





6154 6156 6158 6160 5708 57106154 

Wavelength (A) 



Si 6.3 
Si 6.8 
Si 7.3 



6156 



6158 




4480 4482 5527 . 5528 
Wavelength (A) 



5529 



Figure 12. The spectrum near 6156A, a region containing the Si I 6155A 
and the three lines of O I RMT10. Synthetic spectra are shown for the O and 
Si abundances indicated on the figure. 



Figure 13. The Mgl 5528 A and Mgll 4481 A regions showing a com- 
parison of the observed spectrum versus the synthetic spectra for three Mg 
abundances. 



feature give an abundance about 1.5 dex higher abundance, a dif- 
ference attributed to non-LTE effects. 

Na: The ^/-values are taken from the NIST database. Five Nal 
lines were suitable for abundance analysis from the McDonald 
spectrum (Table 3). RMT 4 and 6 give somewhat different results 
but we assign all lines the same weight. 

Mg: The ^/-values for Mgl and Mgll lines are taken from the 
NIST database. Four Mgl lines are listed in Table 3. Figure [T3l 
shows the best-fitting synthetic spectrum for the 5528 A Mg I line. 



The strongest lines not included in Table 3 are the Mg b triplet at 
5167 (blended with Fel), 5172 (blended with Fel), and 5183 A 
have measured equivalent widths of 470 mA ,413 mA , and 496 
mA in the McDonald spectrum, respectively; these yield an (un- 
certain) abundance of 7.12 dex in good agreement, however, with 
lines in Table 3. The strong Mg II 4481 A feature (Figure[7?t is well 
reproduced by the abundance from the Mg I lines and shows the 
asymmetry in the blue wing attributed to a wind. 

Al: The resonance lines at 3944 A and 3961 A are detectable. Their 
^/-values are from the NIST database. Comparison of observed 
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Table 3. Species used in the analysis and corresponding abundances for a model of T e g = 6500 K, log g = 0.5. RMT is the multiplet number given by the 
Revised Multiplet Table (Moore 1945). 



Species 


A 


EW loge(X) 


LEP 




H1V1 1 




\A) 






(eVJ 






G I 


a 'ion c o n 
4700.02U 


a n 
4U 


Q OQ 
O.Z6 


7.48 


—2.62 


6 


G I 


/i77n nnn 
4 ( /U.UUU 


Z { 


i nn 

/ .yu 


7.48 


—2.44 


6 


G I 


/i 771 i" nn 

4771.690 


1 DC 

1U5 


Q 1 £ 


7.49 


— 1.87 


6 


G I 


,i77cr 07n 

4775.870 


51 


8.10 


7.49 


—2.30 


6 


G I 


4932.049 


ft A 
04 


(.to 


7.68 


—1.66 


13 


G I 


5380.340 


bb 


7.85 


7.68 


— 1.62 


11 


G I 


TI 1 1 A TO 

7111.470 


c c 
bb 


7.90 


8.64 


— 1.09 


26 


C I 


7113.180 


DO 


7.80 


8.65 


—0.77 


26 


C i 


71 1 r i 7n 

7115.170 


c c 
oo 


7.77 


8.64 


—0.94 


26 


G i 


7116.990 


bb 


*7 on 
/ .oU 


8.65 


—0.91 


26 


G 1 


71 i n ftftn 


bb 


(.80 


8.64 


— 1.15 


26 


iv 1 


Q01 ft Q/1 n 

5Z10.O4U 


bb 


/ 7 1Q 

i . iy 


i n qa 


(115 

U.13 


2 


AT t 
JV I 


5ZZo. loU 


bb 


^ /.iy 


10.33 


—0.27 


2 


[Oi] 


tr. CT 77 Q /I n 

007 f .o4U 


C C 
Oj 


O.O^ 


1.97 


—8.23 


3 


[0 1] 


conn ori/i 
ooUU.oU4 


1 Q 


O.O ( 


0.00 


—9.78 


1 


[0 1] 


/;■ o / ■ o 7on 

oooo.7oU 


c c 


o.4Z 


0.02 


— 10.26 


1 


T 

U I 


i?i cr nvn 


bb 


8.60 


10.74 


—0.66 


10 


O I 


6156.760 




8.52 


10.74 


—0.44 


10 


O I 


i^1 ro 1 7A 

6158.170 


bb 


o.o / 


10.74 


—0.30 


10 


iVa I 


A fl70 cr /I 1 

4978.541 


28 


6.15 


2.10 


— 1.22 


9 


Na I 


cr on ijqo 

5682.633 


c c 
bb 


6.02 


2.10 


—0.71 


6 


TV a I 


cr o o on cr 

5688.205 


bb 


t; Q7 
0.5 / 


2.10 


—0.45 


6 


iVa I 


O 1 O O n CT CT 

8183.255 


bb 


D.Z / 


2.10 


0.24 


4 


Na I 


o i n a 7nn 

8194. 79U 


bb 


D.Z I 


2.10 


—0.46 


4 


A/f T 

Mgl 


/I AC7 rflK 


bb 


t .ZZ 


A QC 


—0.90 


16 


Mgl 


a 7no nni 

4702.991 


one; 


7 in 


4.35 


—0.44 


11 


Mgl 


5528.410 


217 


7.26 


4.35 


—0.50 


9 


Mgl 


cr 71 i nnn 

5711.090 


/I V 

47 


6.98 


4.34 


— 1.72 


8 


Mg II 


4481.126 


c o 
DO 


7.17 


8.86 


0.75 


4 


Mg II 


4481.150 


bb 


/ .1 / 


8.86 


—0.55 


4 


Mg II 


4481.325 


c c 
oo 


717 
/ .1 / 


8.86 


0.59 


4 


A 1 i 

z4i I 


on /i /i nnn 

3944.020 


bb 


4.91 


0.00 


—0.64 


1 


A 1 T 

/U I 


on/? i cr no 

3961.523 


c c 


A Q1 

4.ol 


0.01 


—0.34 


1 


5« I 


cr /? cr cr crn 
ODD0.55U 


1 o 
1Z 


*7 nft 


4.92 


— 1.94 


10 


OI I 


r7ni 1 no 
o rUl.lUo 


1 o 


7 on 
i.zy 


4.93 


— 1.95 


10 


bi I 


cr 7no ait 

5708.437 


/I /? 

4b 


7 17 

(At 


4.95 


— 1.37 


10 


bl I 


CT 770 1 Cn 

5772.150 


on 
ZU 


7 1/1 
/ . 14 


5.08 


— 1.65 


17 


C' T 

Dl I 


tT7no n7n 

5793.070 


29 


7.50 


4.93 


— 1.96 


9 


Si I 


6145.020 


c c 


^ 6.99 


5.62 


—0.84 


29 


C" Y 

OI I 


6155.140 


66 


6.90 


5.62 


—0.30 


29 


Si II 


rnrr n o n 

5055.980 


c c 
bb 


6.70 


10.07 


0.52 


5 


C" YY 

bi II 


i" o a 7 nn "i 

6347.091 


285 


8.03 


8.12 


0.15 


2 


Sin 


c o 71 o c n 

6371.359 


252 


7.91 


8.12 


—0.08 


2 


C T 

b I 


/i /^n a i 1 o 

4694.1 lo 


A O 

■Li 


6.98 


6.53 


— 1.71 


2 


O Y 

o I 


a cc\ cr /i /i o 

4o95.44o 


39 


7.08 


6.53 


— 1.87 


2 


S l 


4696.252 


16 


6.88 


6.53 


—2.10 


2 


O Y 

1 


n fT O COO 

6052.583 


c c 
bb 


7.00 


7.86 


— 1.26 


10 


S l 


6743.580 


25 


7.20 


7.87 


— 1.07 


8 


Si 


6748.840 


SS 


6.90 


7.87 


-0.64 


8 


Si 


6757.170 


ss 


6.83 


7.87 


-0.35 


8 


Gal 


4425.441 


95 


5.36 


1.88 


-0.36 


4 


Cal 


4434.957 


SS 


5.50 


1.88 


-0.01 


4 


Cal 


4435.679 


ss 


5.50 


1.88 


-0.52 


4 


Cal 


4578.550 


19 


5.21 


2.52 


-0.56 


23 


Cal 


5581.965 


55 


5.60 


2.52 


-0.71 


21 


Cal 


5588.749 


98 


5.35 


2.53 


0.21 


21 


Cal 


5590.114 


55 


5.60 


2.52 


-0.71 


21 


Cal 


5594.462 


80 


5.45 


2.52 


-0.05 


21 


Cal 


5601.277 


33 


5.59 


2.53 


-0.69 


21 


Cal 


5857.454 


73 


5.47 


2.93 


0.23 


47 


Cal 


6122.219 


131 


5.57 


1.89 


-0.32 


3 


Cal 


6162.172 


161 


5.60 


1.90 


-0.09 


3 


Cal 


6462.566 


122 


5.44 


2.52 


0.31 


18 



Species 


A 


EW loge(X) 


LEP 


l °g{gf) 


H1V1 1 




( A ) 




l^ClCA J 


(eVJ 






Call 


cr n 1 n n7i 

5019.971 


bb 


c; /in 
0.4U 


7.52 


—0.28 


15 


Call 


cr no cr n/" 1 

5285.266 


14 


5.43 


7.51 


— 1.18 


14 


Call 


cr on7 nn ^ 

5307.224 


c c 


K crn 
O.OU 


7.52 


—0.90 


14 


Ca II 


cr o on "i nn 

5339.190 


12 


5.25 


8.44 


—0.33 


20 


C ~ YY 

AcII 


cr no i m n 

5031.019 


i an 

ioy 


o nfi 
z.Uo 


1.36 


—0.40 


23 


C ~ YY 

bC II 


cr non ono 

5239.823 


7Q 


1 QO 


1.46 


—0.77 


26 


bC II 


cr o 1 o o crn 

5318.350 


19 


2.23 


1.36 


—2.01 


22 


Sell 


cr cr n/" 1 o i n 

5526.810 


1 d A 
104 


i nc; 

i.yo 


1.77 


0.02 


31 


Sen 


cr n a n n7i 

OD4U.971 


fin 

oy 


z.lo 


1.50 


— 1.13 


29 


OCII 


cTftftv 1 /in 
ODD / .149 


bb 


1 f\A 


1.50 


— 1.31 


29 


C TT 
OCII 


cr /qt?n n /i n 
0DOy.U4U 


P.A 

04 


O 1 cr 

z.10 


1.50 


— 1.20 


29 


Sen 


cr n o /i "i nn 

0084.190 


nn 

yu 


O Of^ 
Z.ZO 


1.51 


— 1.07 


29 


I % I 


zinnn rrno 

4yyy.oUo 


QQ 

OO 


o.y ( 


0.83 


0.25 


38 


I % I 


oUoo.4oo 


bb 




1 A A 

1.44 


0.13 


110 


111 


cr i no ot;n 

oiyi.yoy 


bb 


^; 4.10 


0.02 


— 1.01 


A 

4 


Ti II 


(ini cr i o 

4493.513 


70 


4.10 


1.08 


—3.02 


18 


Ti II 


a cr /i a nnn 

4544. 009 


i no 
LUZ 


/i n/i 
4.U4 


1.24 


—2.58 


60 


Ti II 


4568.312 


78 


4.19 


1.22 


—2.94 


60 


Ti II 


A 7nO / ■ /"> O 

4708.663 


1 KQ 

lOo 


4. 10 


1.24 


—2.34 


48 


Ti II 


,1 7no cr o cr 

4798.535 


120 


4.11 


1.08 


—2.68 


17 


Ti II 


a o c cr o n 

4ooo.o20 


1 OQ 
LZo 


A QA 

4.oU 


1.11 


—2.79 


29 


Ti II 


/I O 7/1 m A 

4o / 4.U14 


11/1 
1 14 


q nfi 


3.09 


—0.80 


114 


Ti II 


/i n 1 i 1 no 

4yii. iyo 


1 /I Q 

14o 


a riQ 


3.12 


— 0.61 


11/1 
114 


Ti II 


cr n7n no 1 

5072.281 


1 1 9 

llo 


A OA 


3.12 


1.06 


113 


Ti II 


cr i o cr nnn 

5185.900 


173 


4.03 


1.89 


— 1.49 


86 


Ti II 


cr o o i m cr 

5381.015 


bb 


3.73 


1.57 


— 1.92 


69 


Ti II 


r (i i o 7 r" "i 

5418.751 


92 


3.63 


1.58 


—2.00 


69 


V II 


^no/^ 77n 

4036.779 


c c 
bb 


3.14 


1.48 


— 1.53 


9 


V ii 


. i non c 7 a 

4039.574 


c c 
bb 


^ 3.14 


1.82 


—2.11 


32 


V ii 


A CT£? /I CT 70 

4564.578 


bb 


Q Q/1 

o.o4 


2.27 


— 1.39 


56 


Or I 


/i cr a cr n /i cr 

4545.945 


25 


4.80 


0.94 


— 1.38 


10 


Gr I 


,1 C A C 1 A A 

4o4o.l44 


01 


A C7 

4.0 / 


1.03 


—0.71 


21 


Gr I 


a n cr n 1 tro 

4652.158 


29 


A cr 7 
4.57 


1.00 


— 1.03 


21 


Gr I 


cr nn /i cr i o 

5204.518 


1 A O 

14z 


/I 71 
4. / 1 


0.94 


—0.21 


7 


Gr I 


cr nno a nn 

520o.4z0 


197 


4.82 


0.94 


0.16 


7 


Gr I 


cr o c\ f r-o c 

5296.686 


c c 

DO 


/I OI 

4.51 


0.98 


— 1.41 


18 


Cr I 


cnn7 • y / 1 / \ 

5297.360 


c c 

DO 


4.81 


2.90 


0.17 


94 


Cr I 


crnno n c c\ 

5298.269 


DO 


4.71 


0.98 


— 1.16 


18 


Gr I 


C O A O O 1 n 

5348.312 


18 


4.56 


1.00 


— 1.29 


18 


Gr I 


cr /i nn 7nn 

5409.790 


68 


4.71 


1.03 


—0.72 


18 


Cr II 


J o "i n o c n 

4812.350 


68 


4.42 


3.86 


— 1.80 


30 


Gr II 


/I O O A £C\'~r 

4884.607 


CT7 

/ 


4.58 


3.86 


—2.08 


30 


Gr II 


5246.768 


35 


4.56 


3.71 


—2.46 


23 


Gr II 


cr on cr o crn 

5305.850 


c;7 
/ 


/I cr. cr 
4.00 


3.83 


—2.08 


24 


Gr II 


cr ono a a n 

5308.440 


53 


/I /I /I 
4.44 


4.07 


— 1.81 


43 


Gr II 


cr o "i n c o 7 

5310.687 


28 


4.57 


4.07 


—2.27 


43 


Gr II 


c o 1 o c nn 

5313.590 


98 


4.69 


4.07 


— 1.65 


43 


Gr II 


cr /i 70 o/ -1 cr 

5478.365 


52 


4.62 


4.17 


— 1.91 


50 


Gr II 


5510.702 


27 


4.51 


3.82 


-2.45 


23 


Mm 


4041.361 


55 


4.48 


2.11 


0.29 


5 


Mm 


4754.040 


43 


4.36 


2.27 


-0.09 


16 


Mm 


4783.420 


64 


4.47 


2.29 


0.04 


16 


Col 


3995.306 


55 


4.26 


0.92 


-0.22 


31 


Col 


4121.318 


55 


^ 4.06 


0.92 


-0.32 


28 


Nil 


4714.408 


78 


5.34 


3.38 


0.23 


98 


Nil 


4829.016 


26 


5.43 


3.54 


-0.33 


131 


Nil 


4937.341 


19 


5.40 


3.61 


-0.40 


114 


Nil 


4980.166 


42 


5.50 


3.61 


-0.11 


112 


Nil 


5035.357 


55 


5.39 


3.64 


0.29 


143 


Nil 


5080.523 


55 


5.39 


3.66 


0.13 


143 


Nil 


5081.107 


55 


5.54 


3.85 


0.30 


194 


Nil 


5082.354 


55 


5.69 


3.66 


-0.54 


130 


Nil 


5084.081 


55 


5.69 


3.68 


0.03 


162 
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Table 4. (Continued.) 



Species 


X 


EW loge(X) LEP 


i°g(s/) 


RMT 


Species 


X 


EW loge(X) LEP 


log(s/) 


RMT 




(A) 


(mA) 


(dex) 


(eV) 








(A) 


(mA) 


(dex) 


(eV) 






Nil 


5099.927 


SS 


5.44 


3.68 


-0.10 


161 


Yu 


5119.111 


SS 


sC 0.81 


0.99 


-1.36 


20 


Nil 


5155.762 


ss 


5.44 


3.90 


-0.09 


210 


Y u 


5123.210 


SS 


0.81 


0.99 


-0.83 


21 


Cul 


5105.537 


SS 


sC 3.35 


1.39 


-1.52 


2 


Yll 


5200.413 


SS 


sC 0.68 


0.99 


-0.57 


20 


Znl 


4680.138 


19 


3.91 


4.00 


-0.86 


2 


Zr II 


4050.329 


SS 


sS 1.24 


0.71 


-1.06 


43 


Znl 


4722.153 


44 


3.88 


4.03 


-0.39 


2 


ZrU 


4208.980 


SS 


1.14 


0.71 


-0.51 


41 


Znl 


4810.534 


53 


3.80 


4.08 


-0.17 


2 


ZrU 


4211.880 


44 


1.35 


0.53 


-1.04 


15 


Sr 11 


4077.714 


SS 


2.54 


0.00 


0.14 


1 


Zr 11 


4496.960 


54 


1.44 


0.71 


-0.89 


40 


Sr 11 


4215.524 


SS 


2.49 


0.00 


-0.18 


1 


Ball 


5853.675 


79 


1.05 


0.60 


-1.01 


2 


Yu 


4854.867 


SS 


0.88 


0.99 


-0.38 


22 


La 11 


5114.560 


SS 


sC 0.31 


0.24 


-1.03 


36 


Y u 


4883.690 


83 


0.59 


1.08 


0.07 


22 


Ndu 


4303.580 


SS 


< 0.44 


0.00 


0.08 


10 


Yu 


5087.420 


SS 


0.58 


1.08 


-0.17 


20 


Eull 


4129.720 


SS 


-0.05 


0.00 


0.22 


1 



Table 5. Abundances of the observed species for IRAS 18095+2704 are presented for a model atmospheres of T e g = 6700 K, log g = 0.3 (KT), £ = 6.0 and 
T ef j = 6500 K, log g = 0.5, £ = 4.7 (this work). The solar abundances from Asplund et al. (2009) is used to convert Klochkova's abundance of element X to 
[X/Fe] for comparison purposes. 



Species 


K95 


K95 


This work 


l°g e 


loge(X) 


[X/Fe] 


logepsT) 


[X/Fe] 


A 


N 


Cl 


8.27 


+0.50 


7.92 


+0.40 


-0.35 


11 


8.43 


Nl 


7.66 


+0.47 


< 7.19 


sC +0.27 


-0.47 ^ 


2 


7.83 


Oi 


8.74 


+0.65 


8.55 


+0.77 


-0.19 


6 


8.69 


Nal 


6.02 


+0.47 


6.12 


+0.79 


+0.10 


5 


6.24 


Mgl 


7.42 


+0.62 


7.14 


+0.45 


-0.28 


4 


7.60 


Mg 11 


7.52 


+0.72 


7.17 


+0.48 


-0.35 


1 


7.60 


All 


5.81 


+0.12 


4.86 


-0.68 


-0.95 


2 


6.45 


Sil 


7.48 


+0.71 


7.18 


+0.58 


-0.30 


6 


7.51 


Sill 






6.70 


+0.10 




1 


7.51 


Si 


6.96 


+0.53 


6.98 


+0.77 


+0.02 


7 


7.12 


Cal 


5.84 


+0.26 


5.48 


+0.05 


-0.36 


13 


6.34 


Call 






5.40 


-0.03 




4 


6.34 


Sen 


2.25 


-0.30 


2.08 


-0.16 


-0.17 


7 


3.15 


Til 


5.68 


+1.44 


3.97 


-0.07 


-1.71 


1 


4.95 


Till 


4.05 


-0.19 


4.04 


+0.00 


-0.01 


12 


4.95 


Vn 


3.45 


+0.23 


3.24 


+0.22 


-0.21 


2 


3.93 


Crl 


6.02 


+1.13 


4.71 


-0.02 


-1.31 


10 


5.64 


Cm 


4.98 


+0.09 


4.54 


-0.19 


-0.44 


9 


5.64 


Mm 


5.12 


+0.51 


4.44 


-0.08 


-0.68 


3 


5.43 


Fel 


6.71 


-0.01 


6.59 


+0.00 


-0.12 


31 


7.50 


Fell 


6.73 


+0.01 


6.59 


+0.00 


-0.14 


11 


7.50 


Col 


5.74 


+1.60 


4.26 


+0.18 


-1.48 


1 


4.99 


Nil 


6.13 


+0.66 


5.48 


+0.17 


-0.65 


11 


6.22 


Cul 


4.06 


+0.63 


< 3.35 


sC +0.07 


-0.71 ^ 


1 


4.19 


Znl 


4.60 


+0.78 


3.86 


+0.21 


-0.74 


3 


4.56 


Sr 11 






2.52 


+0.56 




2 


2.87 


Yll 


1.42 


-0.04 


0.72 


-0.58 


-0.70 


4 


2.21 


Zr II 






1.31 


-0.36 




3 


2.58 


Ball 


1.13 


-0.22 


1.05 


-0.22 


-0.08 


1 


2.18 


La 11 


0.83 


+0.39 


s: 0.31 


sC +0.12 


-0.52 s: 


1 


1.10 


Ndu 


1.56 


+0.84 


s: 0.44 


sC -0.07 


-1.12 ^ 


1 


1.42 


Eull 


0.96 


+1.23 


-0.05 


+0.34 


-1.01 


1 


0.52 



Notes. K95: Klochkova (1995). 

A = log e(X) This work - log e(X) K95 

N is the number of the lines employed in our abundances determination. 
The solar abundances (log e ) are from Asplund et al. (2009). 
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and synthetic spectra is presented in Figure Q5] Excited Al I lines 
were searched for but not surprisingly were undetectable. 

Si: Selection of Si I lines was made starting with the list of lines 
used by Asplund (2000) for his solar abundance determination. As- 
plund's adopted ^/-values come from Garz (1973) with the adjust- 
ment recommended by Becker et al. (1980). Asplund et al. (2009) 
remark that use of a 3D model solar atmosphere and non-LTE cor- 
rections (Shi et al. 2008) do not change his 2000 estimate for the Si 
abundance. 

The (//-values for Si II lines are those recommended by Kelleher & 
Podobedova (2008). Especially prominent in the McDonald spec- 
trum are the lines at 6347 A and 6371 A. These lines provide an 
implausibly high abundance, a value about 1.0 dex higher than that 
from the Si I lines and corresponding to [Si/Fe] ~ +1.4. This abun- 
dance is likely an indication that the lines are not formed in LTE. A 
search for weak Si II lines yielded the line at 5055.98 A from RMT 
5 lines which provides the abundance log e(Si) = 6.7. This is not 
only 1.3 dex less than the value from the 6347 A and 6371 A lines 
but about 0.5 dex less than the abundance from Si I lines. 

S: The (//-values for S I lines were taken from Podobedova et al. 
(2009). Seven lines from three multiplets are easily measurable. 

Ca: The (//-values for the 13 Cal lines in Table 3 are taken from 
the NIST database. New measurements for RMT 3 by Aldenius et 
al. (2009) are smaller by only 0.07 dex, a difference that is ignored 
here. Figure Q3] shows the best-fitting synthetic spectrum for the 
Ca I 4425 A and 4435 A regions. 

The (//-values for three of the four Ca II lines in Table 3 are taken 
from the NIST database with the entry for the 5339 A line from the 
Kurucz database in the absence of an entry in the former database. 
The four lines give consistent results: Figure [T6] shows the best- 
fitting synthetic spectrum for the Ca II lines at 5019, 5285, and 5307 

A. 

Sc: The (//-values for the Sc II lines in Table 3 are exclusively 
from Lawler & Dakin (1989) who combined radiative lifetime and 
branching ratio measurements. 




5020 5284 5286 5308 

Wavelength (A) 



Figure 16. The Call 5019, 5285, and 5307 A regions showing a compari- 
son of the observed spectrum versus the synthetic spectra for three Ca abun- 
dances. 

Ti: The (//-values for Ti II lines are taken from Pickering et al. 
(2001, 2002). The 4999 A Ti I line gives the abundance in Table 3. 
Neutral titanium abundance is constrained with Ti I lines at 5036 
A and 5192 A for which we set abundance limits in Table 3 using 
^/-values from the NIST database. 

V: The (//-values for Vll lines are taken from Biemont et al. 
(1989). The leading lines in the solar spectrum expected in the spec- 
trum of IRAS 18095+2704 are at 4036.77A and 4564.58A . The 
former is detectable but the latter's absence provides an upper limit 
to the V abundance. 

Cr: The (//-values for Cr I lines are taken from the NIST database. 
Sobeck et al.'s (2007) measurements are within ±0.02 dex of NIST 
values for most lines in Table 3 and always within ±0.10 dex. 
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A majority of the selected Cm lines has (//-values in the NIST 
database and for the missing minority we take semi-empirical val- 
ues from the Kurucz linelist. Nilsson et al. (2006) report new mea- 
surements of Cr II (//-values from a combination of radiative life- 
times and branching fractions. The majority of the chosen stellar 
lines in the paper are longward of 4850 A which is the long wave- 
length limit for the lines measured by Nilsson et al. For the seven 
lines in Nilsson et al.'s list between 4000A and 4850A, the mean 
difference between the NIST and their entries for log gf is just 
—0.06. Therefore, we make no adjustment to the NIST (and Ku- 
rucz) entries. 

Mn: The (//-values for Mnl lines are taken from Blackwell- 
Whitehead & Bergemann (2007 when available or otherwise from 
the NIST database. Hyperfine structure was considered for all 
lines with data taken from Kurucfl as discussed by Prochaska & 
McWilliam (2000). 

Fe: The ^/-values for Fe I and Fe II lines are from Fuhr & Wiese 
(2006). Exclusion of three relatively strong Fel lines at 5232 
(EW:224 mA ), 5429 (EW: 220 mA ), and 5446 A (EW:206 mA ) 
changes the Fe abundance only —0.02 dex. Analysis of these lines 
was discussed in Section 4.1. 

Co: The search for Co I lines drew on the tabulation of gf -values 
provided by the NIST database. Few Co I lines are expected to be 
present: the leading candidates are lines at 3995.31, 4121.32, and 
41 18.77 A. The 3995 A line, after allowance for a blending Fe I line 
gives the abundance in Table 3. The other two lines are present but 
blended. 

Ni: The (//-values for the Ni I lines are taken from the NIST 
database. A good selection of Ni I lines is available for an abun- 
dance analysis. 

Cu: Copper through Cul lines of RMT2 is not detectable in the 
spectrum. The strongest line of RMT2 at 5105.54A gives the upper 
limit in Table 3 with the line's (//-value taken from Bielski (1975). 

Zn: Zinc is represented by the three Zn I lines of RMT2 at 4722, 
4680, and 4810 A. The (//-values are from Biemont & Godefroid 

(1980) . 

Sr: The SrII resonance lines at 4077 A and 4215 A are present as 
very strong lines and too strong for a reliable abundance determi- 
nation, i.e., they have EWs of 473 mA and 408 mA , respectively. 

Y: Selection of Y II lines is based on the solar lines judged to be un- 
blended Yll lines in the solar spectrum by Hannaford et al. (1982) 
who provide accurate gf -values. Figure [IT] shows the best-fitting 
synthetic spectrum for the Y II line at 4883 A . 

Zr: Our search for Zr II lines drew on the papers by Hannaford et al. 

(1981) and Ljung et al. (2006) who measured accurate laboratory 
g f -values and conducted an analysis of Zr II lines to determine the 
solar Zr abundance. Figure Q~8] shows synthetic spectra fits to two 
Zr II lines. 

Ba: The Ba abundance is based on the 5854A Ball line, the weak- 
est line of RMT 2. The (//-value adopted is the mean of the 
experimental values from Gallagher (1967) and Davidson et al. 
(1992). Hyperfine and isotopic splittings are taken into account 
from McWilliam (1998). The stronger lines from RMT 2 and the 
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Figure 17. Observed and synthetic spectra for the Y II line at 4883 A . 




4210 



4212 4494 
Wavelength (A) 



4496 



Figure 18. Observed and synthetic spectra for the Zr II lines at 421 1 A and 
4496 A . 



resonance lines (RMT 1) give a roughly 0.7 dex higher Ba abun- 
dance and are not considered further. 

La: The Lall line at 5114 A region was too weak to measure but 
used to set the upper limit: log e(La) ^ 0.31. The gj-value for the 
line is taken from Lawler et al. (2001a). 

Nd: The absent Nd II line at 4303 A is used to set an upper limit: 
log e(Nd) ^ 0.44. The (//-value for the line is taken from Den Har- 
tog et al. (2003). 

Eu: The Eu II resonance lines are 4 129 A and 4205 A were searched 
for in the spectra. The (//-values, hyperfine and isotopic structure 
are taken from Lawler et al. (2001b). Spectrum synthesis of the 
4 129 A line gives the abundance in Table 5. The 4205A line is too 
seriously blended to yield a useful a Eu abundance. 
Mean abundances are summarized in Table 5. Absolute uncertain- 
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Table 6. Sensitivity of the derived abundances to the uncertainties in the 
model atmosphere parameters for the T e g = 6500 K, log g = 0.5 model. 





Alog e 


Species 


AT ^ 


Alog g 








+ 1 jU 


+U. J 


i 1 A 

+ I.U 


1 A 

-l.U 




( K ) 


(cgs) 


^Km s ) 


(Km s ) 


Cl 


+0.03 


+0.05 


—0.05 


+0.07 


Nl 


—0.06 


+0.14 


—0.03 


+0.04 


Oi 


—0.04 


+0.09 


—0.01 


+0.01 


Nal 


+0.09 


—0.07 


—0.01 


+0.01 


Ma T 


+0.10 


—0.07 


—0.02 


+0.03 


Mg II 


—0.04 


+0.16 


—0.11 


+0.18 


Si i 


+0.10 


—0.06 


—0.01 


+0.03 


Sill 


—0.06 


+0.12 


—0.04 


+0.04 


Si 


+0.08 


—0.01 


—0.02 


+0.02 


Cal 


+0.10 


—0.07 


—0.04 


+0.05 


Ca II 


—0.01 


+0.12 


—0.00 


+0.01 


Sell 


+0.08 


+0.13 


A AO 

— 0.U2 


+0.04 


Til 


+0.15 


—0.06 


—0.01 


+0.02 


Ti II 


+0.08 


+0.13 


—0.03 


+0.06 


Vll 


+0.05 


+0.11 


-0.07 


+0.06 


Crl 


+0.14 


-0.06 


-0.02 


+0.03 


Cm 


+0.04 


+0.12 


-0.02 


+0.03 


Mm 


+0.14 


-0.06 


-0.01 


+0.03 


Fel 


+0.11 


-0.08 


-0.03 


+0.02 


Fell 


+0.05 


+0.13 


-0.03 


+0.03 


Col 


+0.15 


-0.06 


-0.09 


+0.13 


Nil 


+0.12 


-0.06 


-0.02 


+0.03 


Cul 


+0.16 


-0.05 


-0.00 


+0.01 


Znl 


+0.13 


-0.05 


-0.01 


+0.03 


Yll 


+0.09 


+0.12 


-0.04 


+0.06 


Zr 11 


+0.09 


+0.13 


-0.02 


+0.02 


Ball 


+0.17 


+0.01 


-0.04 


+0.06 


hall 


+0.13 


+0.10 


-0.00 


+0.01 


Ndu 


+0.15 


+0.06 


-0.02 


+0.02 


Eull 


+0.13 


+0.09 


-0.08 


+0.15 



ties for the abundances arising from uncertainties of the atmo- 
spheric parameters T e g, log g, and £ are summarized in Table 6 
for changes with respect to the model of +150 K, +0.5 cm s -2 , and 
±1.0 km s _1 for representative lines (i.e., the Co entries are based 
on the EW upper limits for the 3995 mA and 4121 mA lines). From 
the uncertainties listed in Table 6, we find the total absolute uncer- 
tainty to be ranging from 0.08 for C I to 0.19 for Mgll by taking 
the square root of the sum of the square of individual errors (for 
each species) associated with uncertainties in temperature, gravity, 
and microturbulent velocity. In light of the line-to-line scatter of 
abundances, the absolute uncertainties, and, more importantly, the 
probability that the star fails to recognize the suite of assumptions 
behind the model atmosphere and the analysis of the lines, speci- 
fication of abundances to ±0.01 dex is surely astrophysical hubris. 
Thus, we cite values to 0.1 dex in discussion of the abundances. 
In the Introduction, we referred to two results of the abundance 
analysis by Klochkova (1995), the only previously reported analy- 
sis for IRAS 18095+2704. The two results referred to are: (i) large 
differences between the abundance from neutral and ionized lines 
of several elements, and (ii) overabundances (relative to iron) of 
lanthanides. Before discussing the star as a representative of post- 
AGB stars, we comment on (i) in light of our analysis. Point (ii) is 
more properly considered in discussing the star as a member of the 
post-AGB class. 

Previous results for abundance differences between neutral and ion- 



ized lines were huge in several cases: i.e., 1.6 for Ti, 1.4 for V, 1.0 
for Cr, and 2.2 for Y but 0.0 dex for Fe due to the fact that it was 
the condition imposed in determining the surface gravity. Mg I and 
Mgll also gave consistent results (—0.1 dex). Our results provide 
much smaller values, i.e., difference in abundances [X/Fe] derived 
from neutral and first-ionized lines of several elements are: —0.1 
(Ti), and 0.2 (Cr) with again 0.0 dex for Fe. We were unable to de- 
tect V I lines. For silicon, the disparity between abundances from 
Si I lines and the single weak line of Sill from RMT 5 is —0.5 
dex. For Ca, it is ±0.1 dex. We attribute elimination of the ear- 
lier large differences between the abundance from neutral and ion- 
ized lines to application of stringent requirements for the identifica- 
tion of weak lines of neutral species. What about non-LTE effects? 
The fact that the strongest lines generally give abundances that are 
too large according to a considerations of astrophysical plausibility 
is likely due to non-LTE effects: examples include the Oi 7770 A 
triplet, the Si II 6347A and 6371 A lines and the resonance lines of 
SrII and Ball. 



6 IRAS 18095+2704: ABUNDANCE ANOMALIES AND 
RELATIVES 

While IRAS 18095+2704 may be a representative post-AGB star, it 
is not a textbook example of a star in this evolutionary state where 
the AGB star experienced thermal pulses that enrich the envelope 
in products of He-shell burning, i.e., carbon directly from the 3a- 
process and heavy elements from the s-process probably powered 
by the 13 C(q, n) 16 reaction. Then, after a series of thermal pulses 
that provide for a C-rich star with enhancements of s-process nu- 
clides, the star experienced severe mass loss and became a post- 
AGB star with effective temperature increasing over a few thousand 
years. AGB stars of intermediate mass may experience H-burning 
at the base of the deep convective envelope ('hot bottom burning') 
such that the C-rich envelope (with s-process enhancements) is re- 
converted to an O-rich but N-rich envelope before mass loss sets it 
on its post-AGB course. In several ways, IRAS 18095+2704 fails 
to fit these scenarios. 

Our assessment of the evolutionary status of IRAS 18095+2704 
involves (i) identifying abundance anomalies with respect to the 
composition of an unevolved or much less evolved star of the same 
iron abundance of [Fe/H] = —0.9 as obtained from samples of 
halo and (thick) disk stars (see, for example, Reddy et al. 2006), 
and (ii) identifying highly evolved and likely also post-AGB stars 
with compositions similar to those of IRAS 18095+2704. 

We begin by considering those elements whose atmospheric abun- 
dance is expected to be unaffected during the course of stellar evo- 
lution. Such elements are represented in Table 3 and Table 4 by the 
run from Mg to Zn, although Mg and Al are possibly affected by 
internal nucleosynthesis. In the Mg to Zn sequence, the so-called 
Q-elements are expected to be overabundant relative to Fe by com- 
parable amounts, say [a/Fe] ~ ±0.3, provided that the star is truly 
metal-poor. For IRAS 18095+2704, the a-elements have [a/Fe] = 
+0.5 (Mg), +0.6 (Si), +0.8 (S), 0.0 (Ca), and 0.0 (Ti). At first glance, 
these values seem at odds with the roughly equal ratios for these el- 
ements according to studies of disk and halo stars; Mg, Si, and S are 
more enhanced and Ca and Ti are less enhanced relative to Fe than 
in the reference stars with [Fe/H] ~ —0.9. One should not over- 
interpret this result by overlooking potential systematic errors. It 
may be significant that [a/Fe] is roughly correlated with the atom's 
ionization potential and, if this is the case, it may be a signature 
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of systematic errors in either our analysis and/or in the analyses 
of the very different stars including main sequence stars providing 
the reference [a/Fe] values. Across the Fe-group, say Sc to Zn, the 
ratios [X/Fe] for IRAS 18095+2704 are normal to within the ob- 
servational uncertainties, that is [X/Fe] falls within the range ±0.2 
dex. Just possibly, Mn appears overabundant at [Mn/Fe]= —0.1 
where a value of —0.4 is found for disk and halo stars of [Fe/H] 
= -0.9 (Reddy et al. 2006). 

Among the suite of elements thought to be immune to internal nu- 
cleosynthesis is Al for which we find [Al/Fe] to be -0.7. This is 
far from the expected value of about 0.0. Perhaps, the degree of 
ionization of Al departs from LTE values and Al atoms are over- 
ionized relative to LTE. In some post-AGB stars, elements like Al 
which condense readily into grains are underabundant in the stel- 
lar atmosphere (see, for example, Giridhar et al. 2005), and, then, 
the elemental underabundances correlate fairly well with the con- 
densation temperature. This is not the case for IRAS 18095+2704. 
For example, Zn with its 726 K condensation temperature (Lodders 
2003) is underabundant by [Zn/H] ~ -0.7 but S at 680K and Na 
at 95 8K are barely underabundant, say [S/H] ~ [Na/H] ~ -0.1. In 
a few post-AGB stars, the underabundances correlate with the neu- 
tral atom's ionization potential (Rao & Reddy 2005) rather than the 
condensation temperature but again IRAS 18095+2704 does not fit 
this pattern. 

It remains to discuss light elements C, N, O, and Na and the heavy 
elements Y, Zr, Ba, La, and Nd which are candidates for abundance 
adjustments by nucleosynthesis in the course of stellar evolution. 
Clearly, the atmosphere is O-rich (C/O ~ 0.25 by number). Yet, the 
C abundance suggests some C enrichment because [C/Fe] < is 
expected following the first dredge-up but [C/Fe] = +0.4 is found 
here. The O enhancement ([O/Fe] = +0.8) is greater than reported 
for normal stars but the excess seems to mirror that found for S, 
another element of higher than average ionization potential for an 
a-element and the ratio [O/S] ~ is similar to that for normal stars. 
The upper limit to the N abundance is approximately consistent that 
expected from the first dredge-up. In short, the C, N, and O abun- 
dances suggest mild C-enrichment on the AGB following dilution 
of C and enhancement of N during the first dredge-up as the star 
became a red giant following the main sequence (see, for example, 
Iben & Renzini 1983). Textbook enrichment of C by thermal pulses 
with subsequent C-destruction by hot bottom burning is ruled out 
by the lack of appreciable N enrichment. Sodium stands apart from 
this picture: [Na/Fe] = +0.8 is neither a ratio found among nor- 
mal stars nor easily accounted for on nucleosynthetic grounds. The 
Na abundance appears to be affected by non-LTE effects. Non-LTE 
analysis of Na I lines by Takeda & Takada-Hidai (1994) indicates 
the non-LTE corrections (Aloge(X)f\in Nal lines at 4979, 5683, 
5688, and 8195 A is -0.06, -0.15, -0.20 and -0.94 dex respectively 
at temperature 6000 K and gravity 0.5 dex. 

Operation of the s-process would be expected to lead to overabun- 
dances of Y and Zr but underabundances (relative to Fe) are found. 
At [Fe/H] ~ —1, Y and Zr underabundances are not found among 
normal stars. Thus, the Y and Zr abundances for IRAS 18095+2704 
present a puzzle. It is not clear if the same puzzle is provided by 
the heavier elements Ba, La, and Nd for which the s-process would 
provide some enrichment. Barium is nominally consistent with the 
Y and Zr in suggesting a slight underabundance ([Ba/Fe] = —0.2). 
Our search for La and Nd proved unsuccessful and the upper limits 

4 Mogt(X)=\og£(X) NLTE -loge(X) iT £ 



[La/Fe] ~ [Nd/Fe] ~ are consistent with results for normal stars 
and also with the Y, Zr, and Ba abundances for IRAS 18095+2704. 
The r-process Eu abundance is that expected for a normal star; any 
s-process contribution to Eu is expected to be very slight even had 
the s-process operated in IRAS 18095+2704 (Sneden et al. 2010). 
In summary, Y and Zr underabundances represent a puzzle. One is 
struck by the fact that Y and Zr with Al have condensation tem- 
peratures among the highest of the elements in Table 4. All are un- 
derabundant relative to expectation. We noted above that the abun- 
dances do not always correlate well with condensation temperature. 
But elements with condensation temperatures hotter than 1550K 
generally appear among the most underabundant - the set includes 
Al, Ca, Sc, Ti, Y, Zr, and Ba. (La and Nd also fall in the set but 
upper limits to their abundances restricts their relevance here.) At 
[Fe/H] ^ —0.3, the Ca and Ti abundance should be judged with 
respect to their overabundance relative to Fe by about 0.3 dex in 
normal stars, i.e., add about —0.3 dex to the entries for [X/Fe] in 
Table 5 when judging abundance anomalies according to condensa- 
tion temperature. Then, the [X/Fe] for Al, Ca, Ti, Y, Zr, and Ba are 
quite similar for these elements of a similar condensation tempera- 
ture. Scandium appears to be mildly overabundant with respect to 
the speculation that elements of the highest condensation tempera- 
ture (To > 1500 K) are underabundant in IRAS 18095+2704. 

In spite of the above remarks about abundances and conden- 
sation temperatures, similarities with the compositions of some RV 
Tauri variables are present. Among RV Tauri variables and the pre- 
sumably closely related W Vir variables (Maas et al. 2007), the cor- 
relation between abundance anomalies and condensation tempera- 
ture runs from strong to weak. IRAS 18095+2704 definitely falls 
among the latter group. There is a fair correspondence between the 
composition of IRAS 1805+2704 and the RV Tauri variable Al Sco 
(Giridhar et al. 2005) (Figure [79t PI In part, Figure [T9] may reflect 
the fact that systematic errors (e.g., non-LTE effects) are of similar 
magnitude for the two stars with similar atmospheric parameters: 
(T off , log g, [Fe/7J])=(6500, +0.5, -0.9) for IRAS 18095+2704 
and (5300, 0.25, -0.7) for Al Sco. Stars exhibiting a striking cor- 
relation involving condensation temperature seem certain to be bi- 
naries with a circumbinary disk providing infall of gas but not dust 
onto the star responsible for the anomalies (Van Winckel 2003). 
For those stars (e.g., Al Sco) with hints of a correlation involv- 
ing the condensation temperature, their binary status is unknown. 
Certainly, the lack of a strong radial velocity variation for IRAS 
18095+2704 may suggest that it is a single star. But, perhaps, a 
wind off the star, as may be suggested by the strong blue asym- 
metry for strong lines, provides the site for dust-gas separation. In 
summary, we suppose that IRAS 18095+2704 may be related to a 
RV Tauri variable, probably one that has evolved to hotter temper- 
atures beyond the instability strip. 



7 CONCLUDING REMARKS 

Abundance analysis of IRAS 18095+2704 classified as a proto- 
planetary nebula by Hrivnak et al. (1988) suggests the star left the 
AGB before thermal pulses had the opportunity to enrich the at- 
mosphere in the principal products from He-shell burning (C and 
s-process nuclides). The star is not exceptional in this regard - see, 

5 There is little resemblance to the composition of EQ Cas and CE Vir 
for which abundance anomalies correlate well with the atom's ionization 
potential. 
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for example, the abundance analyses of post-AGB stars reviewed 
by Van Winckel (2003). If the star's Fe abundance is adopted as a 
reference, the composition of IRAS 18095+2704 is normal except 
for an underabundance of Al, Y and Zr. A speculation was offered 
that elements including Al, Y, and Zr having a condensation tem- 
perature hotter than 1500 K are underabundant. This suggests that 
the atmosphere is depleted in those elements that condense most 
readily into dust grains. Perhaps, the wind which is suggested by 
the pronounced blue asymmetry to the lines of strong lines removes 
grains selectively. 

Although a variety of spectroscopic indicators provide a con- 
sistent set of atmospheric parameters, the high-resolution optical 
spectra offer evidence that IRAS 18095+2704's atmosphere is only 
an approximation to the classical atmosphere and standard LTE 
analysis techniques used to derive the abundances. The line profiles 
require macroturbulent velocities that are supersonic. Strong lines 
suggest by the presence of blue asymmetries that a wind is driven 
off the star. Refinement of the abundance analysis may primarily 
depend on developing an improved understanding of the physics of 
these dilute extended stellar photospheres. 
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